Fluctuation-dissipation relations outside the linear response regime in a 
two-dimensional driven lattice gas along the direction transverse to the driving force 
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We performed numerical experiments on a two-dimensional driven lattice gas, which constitutes a 
simple stochastic nonequilibrium many-body model. In this model, focusing on the behavior along 
the direction transverse to the external driving force, we numerically measure transport coefficients 
and dynamical fluctuations outside the linear response regime far from equilibrium. Using these 
quantities, we find the validity of the Einstein relation, the Green-Kubo relation and the fluctuation- 
response relation. 

PACS numbers: 05.70.Ln,05.40.-a 
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In order to construct a theoretical framework of 
nonequilibrium statistical mechanics, we seek new uni- 
versal relations characterizing nonequilibrium steady 
states (NESS) far from equilibrium. For this purpose, 
we carry out numerical computations through which we 
study the validity of certain fluctuation-dissipation rela- 
tions. In this paper, we consider one such model, a sim- 
ple stochastic nonequilibrium many-body model which is 
called a 'driven lattice gas (DLG)' 0,0. 

Model: Let r\i be an occupation variable defined on 
each site i = (i x , i v ) of a two-dimensional square lattice 
{(ix,iy)\0 < ix < L,0 < i y < L}. The variable t}% is 
1 if the zth site is occupied by a particle, and if it is 
unoccupied. Periodic boundary conditions are imposed 
by setting rji — rjj, where j — (L,i y ) in the case i x = 0, 
and rji = r)j, where j = (i x ,L) in the case i y — 0. The 
array of all occupation variables, {rji}, is denoted r] and 
called the "configuration" . 

The time evolution of r] is described by the follow- 
ing rule: At each time step, randomly choose a nearest- 
neighbor pair and exchange the values of rji and 
rij with the probability c(i,j;rf) = {1 + cxp[/3Q(r] — > 
J7 lJ )]} -1 , where rfi is the configuration obtained from 
tj through this exchange, and /3 = 1/T is the inverse 
temperature with the Boltzmann constant set to unity. 
Q(rj — > rfi) represents the heat absorbed from the heat 
bath as a result of the configuration change r] — > rf^ . The 
total particle number, N — J2i Vi? is conserved through- 
out the time evolution. The density p — N/L 2 is a pa- 
rameter of the model. Hereafter, we regard the unit of 
time to be L 2 time steps, which is the number of time 
steps for which an arbitrary site is chosen once on aver- 
age. We refer to this time as 1 MCS (Monte Carlo step 
per site). 

In this paper, we study a two-dimensional DLG with 



Q(rj -> rj') = H (v') ~ H (rj) - Ej p (rj -> r,' 



(1) 



where E is an external driving force, and Hq(t]) de- 
scribes an interaction between particles, written Ho(r]) = 
— S(i j) ViVji where denotes a nearest-neighbor 

pair. The quantity j p (r) — > tj') is the spatially-averaged 
current, that is, the net number of particles flowing in 



the x direction: j p (rj —* r\' 

%+(i,o)) - v'tii - %)%+(i,o)(i - ^+(1,0))]- In this stud y> 

we fix j3 — 0.5 in order for the system to be far from the 
critical region (Note that the critical temperature of the 
model with E = is j3 c = 1.76.), and choose large values 
of E in order for the system to be far from equilibrium. 

Our aim: In DLGs, regarded as one of the simplest 
classes of nonequilibrium models, statistical properties 
of NESS have been investigated from various points of 
view. Among them, there is an interesting report that a 
large deviation functional of density fluctuations is shape 
dependent Q. In a two-dimensional DLG, although the 
general properties of fluctuations are quite different from 
those of equilibrium states, it was found numerically that 
a fluctuation relation holds where we consider only 
properties along the direction transverse to the external 
force E. This fluctuation relation is a relation among 
density fluctuations, the chemical potential 0, and 
the temperature of the environment. 

In equilibrium cases, the fluctuation relation is closely 
related to fluctuation-dissipation relations, which relate 
dynamical properties of equilibrium fluctuations with 
transport properties in the linear response regime 0. 
Then noting that the fluctuation relation holds in the 
DLG even far from equilibrium Q, we wish to also in- 
vestigate the validity of fluctuation-dissipation relations 
far from equilibrium, and determine if their equilibrium 
forms hold here as well. 

In order to obtain such relations, we directly measure 
transport coefficients and dynamical fluctuations in the 
direction transverse to the driving force E in the two- 
dimensional DLG investigated above. In spite of the 
fact that these measured values differ from those for 
the equilibrium state(i? = 0), we numerically find that 
three fluctuation-dissipation relations, the Einstein rela- 
tion, the fluctuation-response relation, and the Green- 
Kubo relation, seem to be valid even for NESS far from 
equilibrium. 

Einstein relation: In the linear response regime near 
E = 0, the Einstein relation for interacting many-body 
systems is written 



D X = aT, 



(2) 
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FIG. 1: The quantity ln[— (p(t)) E /A] plotted as a function 
of t in the case (p,E) = (0.5, 10), with A = -0.2, L = 32. 
From the slope of the line, — 0.004t + 1.44, D is estimated as 
0.004 = D(2ty/ L) 2 . The inset displays the L dependence of 
D. The triangles and circles correspond to (p, E) = (0.5, 10) 
and (0.5,0), respectively. 



where D is the density diffusion constant, x is the in- 
tensity of density fluctuations, and a is the conductivity 

a 

In order to investigate the validity of (0) in the di- 
rection transverse to the external driving force E (the 
^-direction), we first need to define the density diffusion 
constant D as the coefficient of the diffusion term in the 
evolution equation describing the averaged behavior of a 
density field in this direction. In this paper, we consider 
the case in which we have prepared as the initial state, a 
steady state under the perturbation potential 



V(i y ) = Asm— ^ 



(3) 



which is obtained by adding ^ i]iV(i y ) to H (rj). Then, 
we remove V(i y ) at t = in order to measure the re- 
laxation of the density field p(t). The function p(t) 
is the Fourier transform of the coarse-grained den- 
sity p(i y ). These quantities are defined as p(t) = 
Ef H =iP(iy)sin^, and p(i y ) = ± £ iG n Iy J]u where 

Mly = {(W^l 1 ^ ^ 1 L.ly}. 

In Fig. ^ choosing A as a sufficiently small value, 
ln(— (p(t)) E /A) is plotted as a function of t in the case 
(p,E) = (0.5,10) with L = 32 and A = -0.2. Here, 
( )X, represents the statistical average under the relax- 
ation process. Because exponentially decaying behavior 
of (p{t))^ is observed, D can be estimated from the form 



{p(t)) V E = const.c- 1 ^ 



(4) 



In the inset of Fig. ^ D is plotted as a function of the sys- 
tem size L in the cases E = and E = 10, with p = 0.5. 
Because both values of D seem to converge, we conclude 



FIG. 2: <tT as a function of D\. The triangle, 
square, star, plus, and circle correspond to (p, E) = 
(0.5, 10), (0.4, 10), (0.3, 10), (0.5, 3) and (0.5,0). The error 
bars represent the statistical errors arising from the fitting 
when D is measured. The thin-dotted line represents D\ = 
ctT. The inset displays the L dependence of a. The trian- 
gles and circles correspond to (p,E) — (0.5,10) and (0.5,0), 
respectively. 



that the size L = 40 can be regarded as sufficiently large 
to study the statistical properties of macroscopic quanti- 
ties in our model. It is important to note here that the 
values of D in the case E = 10 arc different from those 
in the case E = 0. 

Next, we define the conductivity a by adding a suf- 
ficiently small perturbativeg driving force e in the y- 
direction. This is realized by adding the term ejt(rf — > rj') 
to Q(r] — > rj') in 0, where 

i 

- Vii 1 - V^Vi+ioM 1 - Vi+iOA))]- ( 5 ) 

Note that in the x-direction, the particles are still driven 
by E. Then, the averaged current J e in the y-direction 
is defined as 



Je=i(i t (T7-7,'))f' e . 

Using this J c , the conductivity a is written 
a = lim — . 



(6) 



(7) 



In the inset of Fig. [21 a is plotted as a function of the 
system size, L, in the cases E = and E — 10 with 
p = 0.5. Note that the values of a in the case E = 
are smaller than those in the case E = 10, and that 
qualitatively this difference is not the same as that seen 
for D. 

We previously measured the intensity of density fluc- 
tuations x = L£((pf) s - {(pe)f) 2 ), where p e = 
Eien e and n e = {(i x ,i v )\l <i x < L,L/2-£/2- 
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FIG. 3: R(t) plotted as a function of C(t) in the case (p, E) = 
(0.5, 10), with L — 32 and A = -0.2 over the interval < 
t < 800 MCS. The line represents R(t) = 0C(t), where T = 
2(f3 — 0.5). The error bars represent these of R(t), because 
the uncertainty on C(t) is smaller. In the inset, In C(t) is 
plotted. The line there represents — 0.00403t — 1.31. 



1 < i y < L/2 + £/2}. (See Fig. 4 in Ref. Q.) Note 
that £ is chosen so that it satisfies £ <C £ <C L where 
£ is a correlation length. Using these values of in 
Fig. |21 we plot aT as a function of D\ in the cases 
(p,J5) = (0.5, 10), (0.4, 10), (0.3, 10), (0.5, 3) and (0.5,0). 
Noting that the thin-dotted line represents D\ = 
we find that even though the values of D, a and x are 
different from those in the equilibrium case, along the 
direction transverse to the driving force E, the Einstein 
relation (0) is valid, within the precision of the numerical 
computations. 

Fluctuation-response relation: We next study in the 
fluctuation-response relation, which is also a represen- 
tative universal relation in the linear response theory. 
Again in this case, we focus on the properties of the sys- 
tem in the direction transverse to the driving force E. 

First, we employ the same procedure as in the measure- 
ment of D to introduce a time dependent response func- 
tion R(t). That is, we prepare the steady state under the 
perturbation V(i y ) and then remove this perturbation at 
t = 0. Because the profile of the coarse grained density 
p(i y ) is changed by the removal of V(i y ), we make this 
change explicit by defining R(t) in the following form: 



(8) 



We remark that the decaying behavior of R(t) in the case 
E = 10 is plotted as that of - {p{t)) V E / A in Fig. □ 

Next, we introduce the time correlation function of 
density fluctuations in the direction transverse to the 
driving force: 

c(t) = mm)f- (9) 

In the inset of Fig. |3| In C(t) is plotted as a function of 
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FIG. 4: B t t plotted as a function of r in the case (p,E) = 
(0.5, 0) (the dotted line) and (p, E) = (0.5, 10) (the straight 
line) with L = 32. The thin line represents 0.050r and the 
thin dotted line represents 0.055r. In the case E = 10, it is 
seen that the line deviates from the thin dotted line at large 
t MCS. 



time in the case E = 10, with p = 0.5 and L = 32. For 
the NESS, C(t), like R(t) decays exponentially in time. 

Here, in the equilibrium case (E = 0), using R(t) and 
C(t), the fluctuation- response relation is given by 



C(t)=TR{t). 



(10) 



In the NESS far from equilibrium studied here, because 
R(t) and C(t) exhibit a similar behavior, and because 
the fluctuation relation, which is essentially the same as 
C(0) = Ti?(0), has previously been found to hold Q|, we 
conjecture that l|10|) is valid. 

To demonstrate its validity explicitly, in Fig. in the 
case E = 10 with p = 0.5 and L — 32, R(t) is plotted as a 
function of C(t) over the interval < t < 800 MCS. It is 
seen that the slope is equal to 1/T, within the precision 
of the numerical computations. 

Green-Kubo relation: Finally, again considering the 
properties along the direction transverse to the external 
driving force (the y direction) , we investigate the validity 
of the Green-Kubo relation for NESS far from equilib- 
rium. 

Using the spatially-averaged current in the direction 
transverse to the external driving force, jt, defined in JSJl, 
we begin by the r dependent current J T , which represents 
the net number of particles that move in the y direction 
during a time of t MCS, 

1 TL2 

r = -^^2Mv(k-i)^v(k)). (ii) 



k=l 



Then, using this expression for J T , the intensity of the 
current fluctuations is defined by 



(12) 
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FIG. 5: crT as a function of B. The triangle, 
square, star, plus, and circle correspond to (p, E) — 
(0.5, 10), (0.4, 10), (0.3, 10), (0.5, 3) and (0.5, 0). The inset dis- 
plays the L dependence of B. The triangles and circles cor- 
respond to (p,E) = (0.5, 10) and (0.5,0), respectively. 



In Fig^J B t t is plotted as a function of r in the cases 
E = and E — 10, with p = 0.5 and L = 32, respec- 
tively. It is seen that in the case E = 10, the line fitted for 
small times (but much larger than the relaxation time of 
the current correlations) deviates slightly for large times, 
while in the case E = 0, B t t and 0.050r are equal within 
the numerical precision for all times. This bending be- 
havior of B t t might reflect the effect of a long time tail 
in this NESS. 

In the case E = 0, defining B as the slope of B t t, the 
Green-Kubo relation M can be written 



B = oT. 



(13) 



In the case E ^ 0, we define B as the slope of B t t 
obtained from the fitting in the early time regime. With 
this definition, the size dependence of B in the cases E = 
and E = 10 with p = 0.5 is plotted in the inset of Fig. 
EJ The difference between the values of B in the cases 
E = and E = 10 is qualitatively the same as that for 



Considering this similarity between a and B, in Fig. 
13 we plot aT as a function of B in the cases (p, E) — 
(0.5, 0), (0.5, 10), (0.5, 3), (0.4, 10), (0.3, 10) with L = 32. 
The Green-Kubo relation lfT^|l is valid for the NESS con- 
sidered here. However, the deviation seen in Fig. [5] is 
somewhat larger than that in Fig. [21 

Summary: In this paper, we have reported the results 
of numerical experiments on the two-dimensional DLG 
focusing on the properties along the direction transverse 
to the external driving force E. We find that the Einstein 
relation, the fluctuation-response relation and the Green- 
Kubo relation hold in the NESS far from equilibrium. 
(Note that the validity of the fluctuation relation for such 
a state was previously demonstrated in Ref. Q.) 

Compared with this validity of relations in the direc- 
tion transverse to the external driving force, we remark 
the properties along the direction parallel to the exter- 
nal driving force. We studied a one-dimensional DLG, 
and found that the phenomena observed along the direc- 
tion parallel to the external driving force seemed to be 
more complicated than those observed along the direc- 
tion transverse to the external driving force [9(- 

We end with some discussion of the detailed balance 
of fluctuations, which has a deep connection with the va- 
lidity of universal relations in the linear response regime 
near equilibrium states. In our DLG, the detailed balance 
condition for c(i,j;rj) does not hold in the case E ^ 0. 
However, the numerical confirmation of the universal re- 
lations presented here suggests the detailed balance of 
macroscopic fluctuations. We point out that, with re- 
gard to this topic, D. Gabielli et al. studied a stochastic 
model for which the detailed balance condition does not 
hold, and derived the Onsager's reciprocity, which is also 
the linear response relations for macroscopic quantities 

M- 

The author acknowledges S. Sasa and H. Tasaki for 
discussions on NESS. This work was supported by a 
JSPS Research Fellowship for Young Scientists (Grant 
No. 1711222). 



[*] Electronic address: hayashi@jiro.cu-tokyo.ac.jp 
[1] B. Schmittman and R. K. P Zia Phase transitions and 
Critical Phenomena Volume 17, Statistical Mechanics of 
Driven Diffusive Systems (Academic Press, 1995); J. 
Marro and R. Dickman Nonequilibrium Phase Tran- 
sitions in Lattice Models (Cambridge University Press, 
1999) 

[2] S. Katz, J. L. Lebowitz, H. Spohn, J. Stat. Phys. 34, 497 
(1984). 

[3] F. J. Alexander and G. L. Eyink, Phys. Rev. E 57, 6229 
(1998). 

[4] K. Hayashi and S. Sasa, Phys. Rev. E 68, 035104(R) 



(2003). 

[5] This chemical potential was defined operationally. The 
theoretical background employed in [ij is supported by 
the result of Ref. Q . 
[6] S. Sasa and H. Tasaki, |cond-mat/041 1052] 
[7] K. Hayashi and S. Sasa, cond-mat/0507719 
[8] R. Kubo, Rep. Prog. Phys. 29, 255 (1966). 
[9] K. Hayashi and S. Sasa, Phys. Rev. E 71, 046143 (2005). 
[10] D. Gabielli, G. Jona-Lasinio, C. Landim, Phys. Rev. Lett. 
77, 1202, (1996) 



